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Abstract 

The  main  purpose  of  this  paper  is  to  study  the  solution  behavior  of  the  transient 
heat  transfer  problem  for  one  dimensional  symmetric  thermoelectric  shape  memory 
alloy  (SMA)  actuators.  It  is  proved  that  for  the  transient  cooling  problem  with  constant 
electric  current  density  of  magnitude  \J\  there  is  a  value  Jq  of  \J\  such  that  when 
I  J|  <  Jo?  temperature  in  SMA  is  always  decreasing  to  its  steady  state,  and  when 
\J\  >  Jo)  temperature  in  SMA  may  not  be  always  decreasing,  which  is  an  important 
property  of  thermoelectric  SMA  actuators  reported  by  Lagoudas  and  Ding  (1995).  A 
lower  bound  of  Jo  is  given.  The  physical  implications  of  main  results  are  also  discussed. 
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1  Introduction 


Shape  Memory  Alloys  (SMA)  [5]  undergo  an  austenitic  to  martensitic  phase  transfor¬ 
mation  upon  cooling,  while  with  some  hysteresis  they  fully  recover  their  parent  phase  when 
heated  above  the  austenitic  finish  temperature.  A  thermodynamic  formulation  of  the  phase 
transformation  in  SMA  has  been  given  by  Tanaka  and  Nagaki  [14]  based  on  the  thermo¬ 
dynamic  approach  to  dissipation  proposed  by  Edelen  [9].  Many  other  researchers  have  also 
contributed  to  the  thermomechanical  constitutive  description  of  SMA,  some  of  the  key  refer¬ 
ences  can  be  found  in  [2,  3].  SMA  have  been  proposed  as  actuators  for  shape  and  vibration 
control  of  structures,  especially  the  Ni-Ti  (Nitinol)  system.  The  capability  of  producing  large 
actuation  forces  during  the  martensitic  to  austenitic  phase  transition  is  the  major  advantage 
of  SMA.  The  major  disadvantage  in  utilizing  SMA  actuators  is  the  low  rate  of  cooling,  since 
the  time  constant  for  heat  transfer  is  usually  large  compared  to  the  frequencies  required  for 
many  engineering  applications.  Recently  a  new  SMA  thermoelectrically  cooled  or  heated 
actuator  has  been  studied  by  Bhattacharyya  et  al.  [1]  that  utilizes  the  thermoelectric  ef¬ 
fect  for  heat  transfer  from  and  to  the  SMA,  by  making  the  SMA  actuator  the  cold  or  hot 
junction  of  a  thermoelectric  couple  respectively.  From  the  heat  transfer  model  proposed  in 
[1],  Lagoudas  and  Ding  [13]  have  proposed  an  equivalent  simpler  model,  which  is  governed 
by  an  integro-differential  equation,  for  1-D  symmetric  thermoelectric  SMA  actuators.  They 
investigated  both  transient  and  multiple  cycle  solutions  induced  by  a  piecewise  constant 
electric  current  source  numerically.  They  observed  and  conjectured  that  for  the  transient 
cooling  problem  with  constant  electric  current  density  of  magnitude  |  J|  there  is  a  value  Jq 
of  I  J|  such  that  when  |  J|  >  Jo,  the  temperature  in  SMA  may  not  be  always  decreasing. 

The  purpose  of  this  paper  is  to  study  the  behavior  of  transient  solutions  to  heat  transfer 
problems  of  TD  symmetric  thermoelectric  SMA  actuators  theoretically.  We  will  restrict 
ourselves  to  the  linear  case  (the  heat  capacity  of  SMA  is  approximated  by  a  constant)  and 
prove  the  above  important  conjecture  of  the  transient  solution  and  derive  an  estimation  of 
the  value  Jq. 

The  arrangement  of  this  paper  is  as  follows.  In  Section  2,  the  analytical  modeling  is  intro-  . 
duced  and  the  existence  and  uniqueness  of  transient  solutions  are  established.  A  numerical 
example  corresponding  to  the  typical  Nitinol  SMA  actuator  is  also  included  to  display  the 
behavior  of  transient  solutions.  In  Section  3,  we  first  discuss  the  asymptotic  behavior  of 
the  transient  solutions  and  then  prove  that  for  the  transient  cooling  problem  with  constant 
electric  current  density  of  magnitude  jJj  there  is  a  value  Jq  of  |  J|  such  that  when  |  J]  >  Jq, 
the  temperature  in  SMA  may  not  be  always  decreasing.  An  estimation  of  the  value  Jq  is 
also  given. 
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2  Heat  Transfer  Model  of  Thermoelectric  SMA  Ac¬ 
tuators 


P-semiconductor  SMA  N-semiconductor 


X 


Figure  1:  One  dimensional  three  phase  P-SMA-N  thermoelectric  system 


The  1-D  heat  conduction  problem  in  which  a  thin  SMA  material  is  confined  between  two 
semiconductors,  a  N-type  and  a  P-type,  forming  the  junction  of  a  thermoelectric  element 
is  shown  in  Fig.  1.  This  thermoelectric  SMA  element  is  proposed  by  [1]  as  an  efficient 
way  of  fast  cooling  or  heating  of  SMA  actuators.  The  principle  of  thermoelectric  heat 
transfer  is  based  upon  the  Peltier  effect  [8],  according  to  which  an  electric  current  fiow 
creates  a  temperature  differential  between  junctions  where  dissimilar  metals  meet.  For  the 
case  shown  in  Fig.  1,  the  net  heat  flux  across  interfaces,  N-semiconductor/SMA  and  P- 
semiconductor/SMA,  is  given  by 


(Qs  ~  Qn)  ■  ^ 


(2.1) 


(<3s  -  Qp)  • 

where,  for  i  =  N,S  and  P  corresponding  to  A -semiconductor,  SMA  and  P-semiconductor 
phases  respectively,  Qi  is  the  heat  flux  vector  in  the  f-th  phase,  n  is  a  unit  normal  along  the 
positive  X— direction  (as  shown  in  Fig.  1),  a*  is  the  Seebeck  coefficient  of  i— th  phase,  Ti{x,  t)^ 
is  the  temperature  at  {x,t)  in  i-th  phases,  and  J  is  the  current  density.  When  current  J 


flows  through  the  junction,  the  Peltier  effect  induces  jump  discontinuities  in  the  heat  flux 


across  the  interfaces,  which  makes  the  thermoelectric  element  in  Fig.  1  either  a  heating  or 
cooling  element,  depending  on  the  current  direction. 


In  the  1-D  three  phase  N-SMA-P  thermoelectric  system,  it  is  assumed  that  the  tempera¬ 
ture  in  each  phase  does  not  significantly  varies  on  the  plane  perpendicular  to  x— axis,  hence 
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the  Fourier  law  of  heat  conduction  becomes 


dT: 


where  Ki  is  the  thermal  conductivity  of  i-th  phase.  The  electric  current  density  vector 
is  given  by  /(x,t)  =  J{t)n.  If  To  is  the  surrounding  temperature  (for  example,  the  room 
temperature)  of  the  1-D  three  phase  N-SMA-P  thermoelectric  system,  the  convective  heat 
transfer,  which  occurs  across  the  sides  of  the  i— th  phase,  is  included  approximately  as  a 
source  term,  -HP{Ti  -  Tq)/A,  in  the  heat  conduction  equation  for  the  N-SMA-P  system, 
where  H  is  the  heat  convection  coefficient,  P  and  A  are  the  perimeter  and  area  of  the  cross 
section  respectively,  which  are  assumed  to  be  independent  of  x  (see  [6]).  The  heat  conduction 
equations  for  the  three  phases  of  the  N-SMA-P  system  are  then  given  by  (see  [1]) 


Kpf-^^{x,t)  +  —  H—(Tff{x,t)  —To)  —  (x.t),  (2.2) 

d/2  <  X  <  L-fd/2,  t  >  0, 

ffinr  P  &T 

K,^{x,t)+psJHt)-Hj{Z{x,t)-TQ)  =  Ct-^ix,t),  (2.3) 

-d/2  <x<  d/2,  t  >  0, 

Kp^ix,t)  +  peJ^{t)-Hj{Tp(x,t)-To)  =  C^^(x,t),  (2.4) 

-L-dl2<x<  -d/2,  t  >  0. 

where  p,  is  the  electrical  resistivity  and  PiJ^{t)  represents  the  Joule  heat  in  the  i— th  phase, 
and  Ci  is  the  heat  capacity  per  unit  volume  of  i— th  phase.  The  interface  conditions  consist 
of  continuity  of  the  temperature  field 

T,(P)=T„{P),  T.(-|()=7’p(-|0,  (2-5) 

while  the  heat  flux  interface  conditions  (2.1)  given  earlier  reduce  to 

<  ""  (2.6) 
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where  it  is  assumed  that  ot^  =  0,  which  is  justified  as  the  Seebeck  coefficient  of  an  SMA 
is  very  low  compared  to  that  of  the  N  or  P  semiconductor.  Moreover,  the  end  boundary 

conditions  are 


Tn{L  +  -,<)—  To, 

while  the  initial  conditions  are  stated  as 

Tm{x,Q)  =  To, 

^  T,(x,0)  =  To, 

Tp(x,0)  =  To, 


Tp{-L-^,t)=To, 


d  ^  d 


d  d 
"2  -  -  2’ 


d  d 

.L-  -  <x< 

2  -  -  2 


(2.7) 


(2.8) 


The  heat  capacities,  Cf  and  Cf ,  of  N-type  and  P-type  semiconductors  are  taken  to 
be  independent  of  temperature,  while  the  heat  capacity  of  the  SMA  material  varies 
significantly  with  T  due  to  the  temperature  induced  phase  transformation.  Even  though  the 
electrical  resistivity  changes  with  temperature  during  the  phase  transformation  of  SMA, 
lack  of  precise  experimental  data  limits  us  in  assuming  that  /?,  is  constant.  Dependence  of 
Ps  on  temperature  with  possible  hysteresis  can  be  accounted  for  mathematically  by  varying 
the  magnitude  [Jj  of  J. 

Due  to  the  relatively  large  thermal  conductivity  AT,  of  the  SMA,  compared  with  the 
thermal  conductivities  Kp  of  the  F— semiconductor  and  of  the  AT— semiconductor,  the 
temperature  distribution  in  the  SMA  material  is  expected  to  be  almost  constant  for  small 
ratios  djL.  It  is  possible,  therefore,  to  simplify  the  three  phase  model  into  an  approximate 
two  phase  model,  by  incorporating  the  SMA  into  the  P-N  interface.  Integrating  the  field 
equation  (2.3)  with  respect  to  x  from  —d/2  to  d/2,  substituting  (2.6)  in  the  resulting  equa¬ 
tion,  and  then  using  (2.5)  and  the  assumption  that  T,(x,t)  is  independent  of  x,  we  obtain 
the  following  equation 


(2.9) 


Since  the  heat  conduction  equation  (2.3)  for  SMA  is  replaced  by  the  interface  equation  (2.9), 
for  simplicity  of  notation,  we  shift  the  domains  of  the  remaining  field  equations  (2.2)  and 
(2.4)  for  N-type  and  P-type  semiconductors,  [-L-|,  -|]  and  [f,  A  +  f],  to  [-L,0]  and  [0,L], 
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respectively.  If  we  assume  further  that  =  —a^  >  0  and  all  other  material  parameters  of 
P— type  and  N-type  semiconductors  are  the  same,  namely,  =  C„,  Kn  =  Kp  =  K 

and  ps  =  Pp  =  Pi  ^^en  the  temperature  distributions  T^(a;,t)  and  Tp{x,t)  are  symmetric, 
i.e.  Ttf{x,t)  =  Tp{-x,t).  Let  T{x,t)  =  TN{x,t)  =  Tp{-x,t),  the  initial  boundary  value 
problem  for  the  symmetric  1-D  N-SMA-P  system  (2.2)-(2.8)  is  then  reduced  to 


=  K^{x,t)+pJ^{t)  -  ^{T{x,t) -To),  0<x<L,t>0 
T{x,0)=To,  0<x<L,  T(0,0)  =  To, 


<  T{L,t)=To,  t>0,  (2.10) 

2K^{0,t)  +  p,dJ^{t)  -  ^{T{0,t) -To) 

=  -2T{0,t)a^J{t)  +  Ctd^{0,t),  t  >  0. 


Because  of  the  interface  condition  in  equation  (2.10),  classical  methods  such  as  the  sep¬ 
aration  of  variables  method,  integral  transform  methods,  etc.  can  not  be  directly  applied 
to  (2.10)  to  derive  an  analytical  solution.  Furthermore,  because  of  the  temperature  depen¬ 
dence  of  the  heat  capacity  of  SMA,  the  boundary  value  problem  (2.10)  becomes  extremely 
challenging  to  study  theoretically  and  numerically. 

One  way  to  solve  equation  (2.10)  is  expressing  the  interface  equation  in  (2.10)  as  an 
independent  equation  for  T(0,  t),  then  solving  the  following  Dirichlet  boundary  value  problem 
(once  the  distribution  T{0,t)  is  found): 


C„^(x,t)  =  K'^{x,t)+pJ\t)  -  L-{T{x,t) -To),  0<x<L,t>0 
T{x,0)=To,  0<x<L, 

<  (z.lij 

T(L,t)=To,  t>0, 

T(0,f)=T(t),  t>0. 

To  express  the  interface  equation  in  (2.10)  as  an  independent  equation  for  r(0,t),  we  need 

dT 

to  find  the  relation  between  T{0,t)  and  (0,t).  Assume  T{0,t)  =  T{t)  is  given  in  (2.11), 

ox 
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by  using  the  separation  of  variables  method  and  Fourier  series  technique,  w^ould  obtain 
T{x,t)  =  from  (2.11).  Hence  the  relation  between  T(0,t)  and  would 

be  given  by  ^(0,i)  =  Based  on  this  idea,  we  obtain  the  following  integro- 

diflferential  equation  of  T(0,t)  =  T(<),  which  is  equivalent  to  the  interface  equation  in  (2.10) 
(the  detailed  calculation  can  be  found  in  [13]), 


G{t  -  t)  —  Tojr(t) 

r(o)  =  To, 


(2.12) 


where 


and 


CtdL  HP  -La^J{t)  .  1  PHLd 

=  ^i  =  —\  ^2{t)  = - TTT —  +  -  + - 


AK 


C„A’ 


2K 


2  AKA  ’ 


«(,).}  Kfj-ggj’.g).). 


HP  2/5  /•'  X  9/  X  ,  1  PHLd  psdL  .o.^x 

m  =  I  0{t-  r)-^dr  r)J  (r)dr  +  5  +  -^  +  (*) 


Once  the  temperature  distribution  in  the  SMA  interface  T{0,t)  —  T(t)  is  found  from 
(2.12),  the  temperature  distribution T{x,t)  in  P-semiconductor  and  fV-semiconductor  can 
be  found  by  solving  the  initial  boundary  value  problem  (2.11).  Therefore  the  study  of  the 
transient  heat  conduction  for  the  1-D  thermoelectric  P-SMA-N  system  reduces  to  the  analysis, 
of  the  integro-differential  equation  (2.12).  Before  we  proceed  with  further  discussions,  it 
is  desirable,  for  simplicity  of  mathematical  and  physical  discussion,  to  non-dimensionalize 
equation  (2.12)  by  using  non-dimensional  groups  which  are  obtained  by  the  application 
of  the  Buckingham- TT  theorem  [4].  The  non-dimensional  quantities  in  this  study  following 
the  notation  used  in  [1],  where  a  letter  with  an  overbar  represents  the  non-dimensional 


counterpart  of  a  physical  quantity  with  same  letter,  are  listed  below: 


< 


J=JL 


fr  PHL\ 

KA' 


(2.13) 


From  now  on,  unless  mentioned  otherwise,  all  quantities  are  assumed  to  be  non-dimensional. 
For  simplicity  of  notation,  we  will  drop  the  overbar  from  the  letters  used  for  the  non- 
dimensionalized  quantities  introduced  in  (2.13).  Then  equation  (2.12)  becomes 


where 


.  m  =  1, 


(2.14) 


(2.15) 


F(t)  =  y:  « 

n=l 


n=l 


and 


J^{t)  =  H  G{T)dr  +  2  F{t  -  r)J^ir)dr  +  ^  +  ^Hd  -F  ^pdJ\t) 


The  existence  and  uniqueness  of  solution  to  (2.14)  have  to  be  established  first. 

Theorem  2.1  Let  J  6  C(0,oo)  and  the  heat  capacity  of  SMA,  Cl,  be  a  constant.  Then  the 
integro-differential  equation  (2.14)  admits  a  unique  solution  in  C[0,oo). 

Proof.  Let  us  first  consider  the  following  equation 


—{t)  +  HTit)  =  5(0,  t>0 

at 


T{0)  =  1. 


(2.16) 
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We  obtain 


T(()  =  e-^*  +  t  -  ’■)5(T)<iT.  (2.17) 

Substitute  T{t)  into  equation  (2.14),  we  have  the  following  integral  equation, 

fiS{t)  +  G*{t,  T)S(r)dT  =  (2.18) 

Jo 

where 

^"(1)  =  P(t)  -  W)  -  ■ 

Equation  (2.18)  is  a  Volterra  equation  of  the  second  kind,  where  the  kennel  G'{t,T)  has  a 
weak  singularity  of  |t  —  T\~i .  By  [11,  Theorem  7,  pp.  35-37],  equation  (2.18)  has  a  unique 
solution  in  C(0,  oo).  ^ 


Table  1:  Material  parameters 


Ni-Ti  SMA 

P-Element 

Thermal  Conductivity  K 

2.2  X  10“W/(  mm  K) 

1.63  X  10“^W/(  mm  K) 

Heat  Capacity 

2.12X  10-V/(  mm^-°K) 

4.35  X  10-V/(  mm^  K) 

Resistance  Density  p 

6.3242  X  10"^n  -  mm 

1.15  X  10"^n  -  mm 

Seebeck  Coefficient  a 

1.2  X  \Q~^Voltsl°K 

2.15  X  l^-'^Voltsl^K 

Let  us  also  look  at  an  example  of  a  typical  Ni-Ti  SMA  actuator  to  get  some  intuition.  The 
material  parameters  of  the  Ni-Ti  SMA  and  the  P-type  semiconductor  at  room  temperature 
To  =  300°^'  shown  in  Table  1  are  adopted  from  [1].  The  heat  convection  coefficient  H  is 
assumed  to  be  =  2.5  x  10“^W/  mm^  K,  and  the  geometric  parameters  of  the  SMA- 
actuator  are  given  by  L  =  4.0  mm,  d  =  2.0  mm,  P  =  4.0  mm  and  A  =  1.0  mrn?.  By  (2.13), 
we  obtain  the  following  non-dimensionalized  parameters: 

p  =  5.4993  X  10-^  a  =  0.8601,  =  0.4874, 

H  =  0.9816,  d  =  0.25. 

Applying  the  finite  difference  scheme  to  (2.14)  [13],  we  obtain  the  temperature  distribution 
at  the  SMA  interface  corresponding  to  different  values  of  J  for  constant  C^,  which  is  shown 
in  Figure  2. 
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Figure  2:  Temperature  in  SMA  corresponding  to  different  values  of  J  for  constant 

The  temperature  distributions  at  the  SMA  interface  corresponding  to  different  values  of 
J  for  variable  is  shown  in  Fig.  3,  where  the  dimensional  heat  capacity  of  a  typical  Ni-Ti 
SMA  has  been  experimentally  determined  and  may  be  represented  by  the  function  (see  [1, 
12]) 


In  inn  2lal00  | 

Ct  =  C^,  +  [L„-  Cl(M.  -  M,)]  ■  ’  M,<T<M.,  (2.19) 

during  the  forward  transformation  (austenite  to  martesite)  and  the  function 


Ci  =  Cl  +  \L„-Cl{A,-A.)] 


In  100  - 


21nl00  |j.  ^<+^/| 


\At-A • 


Af  <T  <A, 


(2.20) 


during  the  reverse  transformation  (martesite  to  austenite),  with  =  C®  at  all  other  tem¬ 
peratures.  When  an  SMA  is  initially  at  the  austenitic  state,  its  crystal  structure,  upon 
cooling,  starts  to  change  from  the  austenitic  phase  to  the  martesitic  phase  at  the  starting 
temperature,  M^,  and  finishes  at  the  finishing  temperature,  Mf,  when  the  martesitic  phase 
is  formed.  When  an  SMA  is  initially  at  the  martesitic  state,  its  crystal  structure,  upon 
heating,  starts  to  change  from  the  martesitic  phase  to  the  austenitic  phase  at  the  starting 
temperature,  A^,  and  finishes  at  the  finishing  temperature.  A/,  when  the  austenitic  phase 
is  formed.  For  a  typical  Ni-Ti  SM.A.  [1],  M;  =  278'’ A',  M„  =  296° AT,  A^  =  302° A'  and 
Af  =  324° A",  Lii  =  0.1  Jlvnn^  and  C®  =  2.12  x  10"^  J/(  mni^  A'). 
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Figure  3:  Temperature  in  SMA  corresponding  to  different  current  densities  for  variable 

For  the  constant  case,  we  find  from  Fig.  2  that  the  temperature  profiles  in  SMA  are 
quite  different  corresponding  to  different  values  of  the  current  density  J.  For  J  =  — 0.4  and 
J  =  —0.9,  the  corresponding  temperature  in  SMA  always  decreases  asymptotically  to  the 
steady  state  temperature.  For  J  =  -1.5  and  J  =  -1.9,  the  corresponding  temperature  in 
SMA  decreases  first,  it  reaches  its  lowest  value  at  some  finite  time  and  then  asymptotically 
increases  to  the  steady  state  temperature.  Similarly,  for  the  variable  C*  case,  we  observe 
the  same  phenomena  as  shown  in  Fig.  3  and  discussed  in  more  detail  in  [13].  As  pointed 
out  before,  the  major  disadvantage  in  utilizing  SM.4  actuators  is  the  low  rate  of  cooling, 
because  the  time  constant  for  heat  transfer  is  usually  large  compared  to  the  small  time 
constants  required  for  many  high  frequency  engineering  applications.  Thus  to  maintain  a 
monotonic  decay  of  the  SMA  temperature  is  very  important  in  engineering  applications. 

It  is  conjectured  and  confirmed  numerically  in  [13]  that  for  the  transient  cooling  problem 
with  constant  electric  current  density  of  magnitude  |  J|  there  is  a  value  Jq  of  jJj  such  that  - 
when  I  J|  <  Jo,  the  temperature  in  SMA  is  always  decreasing  to  its  steady  state,  and  when 
I  Jj  >  Jq,  the  temperature  in  SMA  may  not  be  always  decreasing.  We  will  further  investigate 
this  question  in  the  next  section. 
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3  Solution  Behavior  of  the  Transient  Heat  Trans¬ 
fer  Problem 


In  this  section,  we  investigate  the  solution  behavior  of  (2.14).  Due  to  the  complicated 
nature  of  SMA  heat  capacity  Q  (see  (2.19)  and  (2.20)),  we  only  study  the  constant 
case.  Thus  through  this  section,  we  always  assume  is  constant.  In  most  of  engineering 
applications  of  thermoelectric  SMA  actuators,  the  typical  time  for  current  density  J{t)  to 
reach  its  steady  state  is  much  smaller  than  the  typical  time  for  temperature  T(t)  to  reach 
its  steady  state.  Therefore  for  the  study  of  long-time  asymptotic  behavior  of  T{t),  we  will 
assume  J{t)  to  be  constant.  For  the  study  of  small-time  asymptotic  behavior  of  T{t),  we  will 
not  assume  J{t)  to  be  constant.  If  J{t)  is  assumed  to  be  constant,  the  exact  solution  of  (2.14) 
can  be  found  by  using  the  Laplace  transform.  However,  the  expression  is  too  complicated 
for  an  immediate  interpretation. 

Let  us  first  look  at  the  long-time  asymptotic  behavior  of  the  transient  solution  to  (2.14). 


Theorem  3.1  Let  J{t)  =  J  be  a  constant  function  in  equation  (2.14)  T(t)  be  the 
solution  of  (2.14).  Then  the  transient  solution  T{t)  is  stable  provided  that 


J  Jmax  — 


-f- 1  -I-  Hd 
a 


The  steady  state  solution  is  given  by  the  following  formula, 

_  {2Ht)i  -b  1  4-  Hd)  -j-  {pd  -f-  47^2) 
{2HT]i  +  l  +  Hd)-aJ 


where  rji  and  i}2  are  given  by 


Vi  =  Jl 


(2n  -  1)^2  H 


(3.1) 


Proof.  To  investigate  the  stability  of  the  solution  to  equation  (2.14),  we  only  need  to  study 
the  transfer  function  of  (2.14),  w{s),  which  is  the  Laplace  transform  of  w{t),  the  solution  of 
(2.14)  corresponding  to  the  unit-impulse  input  at  t  =  0.  Applying  the  Laplace  transform  to 
the  equation 

^  G{t  -  r)  (^(t)  +  Hw(t)^  dr  -|-  M^(0  +  =  ^(0. 

where  6{t)  is  the  unit-impulse  function  at  t  =  0,  we  obtain  the  transfer  function  w{s)  of 
(2.14)  given  by 
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1 


where 


10(5)  = 


G(s)(5  +  H')  +  fJ.S  +  u 


G{s)  =  E 


n=l 


1 

s  +  n^TT^  +  H 


Figure  4:  Graph  of  G{s){s  +  H)  +  fxs  +  u 
For  simplicity  of  notations,  let 


v„  =  n^TT^  +  H,  n  >  1. 

Notice  that  the  denominator  of  w{s),  G{s){s  +  H)  +  fis  +  u,  is  continuous  on  {— Un+i,  —Vn), 
n  =  1,  2  •  •  •,  and  (-ui,  +00);  G(s)(5  +  H)  +  ij,s  +  u  qpoo.  when  5  -*•  u„  ±  0  for  any  n  >  1;’ 
and  ^  ^ 

{G{s){s  +  H)  +  iis  +  u^  =  E  ^  +  fi>0. 

Thus  G(s)(5  +  H)  +  us  +  1/  is  always  increasing.  The  graph  of  G(5)(5  +  H)  +  fxs  -h  u  is 
sketched  in  Fig.  4.  Denote  by  {A„}~^o  zeros  of  G{s){s  +  H)  +  ns  +  uon  the  real  line  3?. 
Thus 


■  •  •  <  A„  <  -Un  <  A„_1  <  -u„_i  <  •  •  •  <  -U2  <  Ai  <  -Vi  <  0, 
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and  A„  ~  0{n^).  It  is  also  easy  to  check  that  G(s)(s  +  if)  +  /xs  +  i/  has  no  other  zero  on 
the  complex  plane.  Thus  except  for  {A„}^o,  w{s)  has  no  other  poles  on  the  complex  plane. 
For  the  first  pole  Aq  of  w{s),  Aq  <  0  if  and  only  if 

6(0)(0  +  H)  +  fi-0  +  u  =  GiO)H  +  u>0. 


By  using  notations  in  (2.15)  and  (3.1),  we  have 

Ht]i  +  ^  +  ^Hd  -  ^aJ  >  0. 

Thus  the  transient  solution  T{t)  of  (2.14)  is  stable  if  and  only  if 

^  2,Htii  “f"  1  d"  //d 

j  ^  ■  • 

a 

When  J  satisfies  the  above  inequality,  the  steady  state  of  the  transient  solution  T{i)  to 
equation  (2.14)  can  be  obtained.  It  is  easy  to  verify  that 

Inn  JF(t)  =  (^T}i  +  ^  +  2"^^)  (2^^ 


dT 

Thus  letting  — (t)  =  0  and  t  00  in  (2.14),  we  obtain  that 
dt 

_  {2Ht]i  +  1  +  Hd)  +  {pd  +  im) 
°°  (2Ht]i  +  1  H"  Hd')  —  ocJ 


I 


Remark  3.1  Note  that  in  the  present  work,  the  heat  convection  coefficient  H  is  assumed 
to  be  constant,  which  is  true  for  small  temperature  variations.  For  the  example  given  in  the 
last  section,  for  a  stable  transient  solution  to  exist,  J  has  to  be  less  than  Jmax  =  1.8056. 
Since  Too  approaches  infinity  when  J  is  close  to  Jmax  or  goes  to  —00,  as  shown  in  Fig.  5,  the 
assumption  of  H  being  constant  is  not  valid  any  longer. 


Remark  3.2  Theorem  3.1  gives  an  important  condition  on  the  current  density  J  for  the 
stable  temperature  of  (2.14).  For  the  transient  cooling  problem  of  SMA  actuator,  i.e.  J  <  0, 
the  transient  solution  is  always  stable.  Corresponding  to  the  example  given  in  the  last  section 
with  the  data  given  in  Table  1,  we  have 


TM  = 


1.5530  +  0.4765J2 
1.5530-  0.8601  J’ 


From  the  graph  of  Too{J),  shown  in  Fig.  5,  we  have  that  =  mmToo(/)  =  0.8285  when 
=  -0.7476. 


14 


(T) 


Figure  5:  The  graph  of  Too 

Even  though  T^'"  is  a  measure  of  the  maximum  ability  of  the  thermoelectric  element  to 
cool  down  the  N-SMA-P  junction,  the  short-time  response  is  also  important  for  the  design 
of  high  frequency  SMA  actuators,  and  it  will  be  investigated  next.  Let  to  >  0  be  very  small, 
and  let  J(t)  be  given  by 


J(t)  —  Jl  o(t),  0  ^  t  <  fo- 

Assume  that  T(t),  for  small  time  t,  has  the  asymptotic  series  expansion  given  by 

T(t)  =  1-1-  T\t  -|-  o(t),  0  ^  t  <  to- 

We  will  determine  T\  from  (2.14)  in  terms  of  Ji,  J2  and  material  parameters.  Substituting 
the  series  expansions  of  T{t)  and  J(t)  into  (2.14),  we  have 

/'  G(t  -t){Ti-VH  +  0(1))  dr  -b  -f  0(1) 

Jo 

+  +  J2t  +  o(t))  +  -  -f  <^(0)  —  -^(O- 


By  comparing  the  coefficients  of  0-th  order  terms  in  t  and  using  the  fact  that 

1 


0  <  F(t)  <  G(t)  < 


t  >  0, 


(3.2) 


we  obtain 
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Ti  = 

where  the  approximate  value  is 
example  given  earlier. 

If  Ji  =  0,  hence  Ti  =  0,  the  above  expansion  of  T{t)  does  not  provide  any  further 
information.  Let  J{t)  be  given  by 

J{t)  =  Jit  +  o(t),  0  <  t  <  to- 

Assume  that  T(t),  for  small  time  t,  has  the  asymptotic  series  expansion  given  by 

T{t)  =  1  +  Tit^  +  o{f),  0  <  t  <  to- 

Substituting  the  series  expansions  of  T(t)  and  J{t)  into  (2.14),  we  have 

G{t  -  r)  (if  +  2T2t  +  o{t))  dr  +  2fiT2t  +  o(f) 

Jo 

+  (-^a{J2t  +  o{t))  +  i  +  i/id)  (1  +  +  o(t2))  =  jr(t). 

By  comparing  the  coefficients  of  first  order  terms  in  t  and  using  (3.2),  we  obtain 

^  ^  r  ^  T 

^2  -  ^-^2  - 


'a  pd  j  aJi  _  qJi 

-  +  -Jt  c.j’ 


justified  whenever 


pdJi 


a 


«  1,  which  is  the  case  for  the 


Next  we  will  study  the  solution  behavior  of  equation  (2.14).  The  following  important 
theorem  explains  why  the  solution  behavior  of  equation  (2.14)  is  quite  different  for  different 
values  of  J. 

Theorem  3.2  LetT{t)  be  the  solution  of  equation  (2.14)  —  J  corresponding 

to  the  transient  cooling  problem.  Then  there  exists  a  value  Jq  for  \J\,  such  that  when  0  < 

1  J]  <  Jg,  T{t)  is  always  decreasing  to  its  steady  state  Too,  otherwise  T{t)  may  not  be  always 
deceasing.  A  lower  bound  of  Jq  is  given  by 

^  —b  +  ^/W+Tac 
=  2^  ’ 

where 
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a  = 


2^/3a 


TT^ 


b  =  pd+^  +  ^[l  +  Hd{l-Ct)], 

TT  TT'^ 


c  =  a. 


Proof.  Let  us  consider  the  homogeneous  equation 


(3.3) 


[  Ti(0)  =  1. 

Applying  the  Laplace  transform  to  this  equation,  we  have 

pL  +  G{s) 


Ti(s)  = 


G(5)(5  +  H)  +  fJLS  +  1/ 

By  Lemma  A,2  in  the  Appendix,  there  exists  a  sequence  of  real  numbers  {6^}  such  that 

^  b 

Ti{^)  =  E7^ 


Then 


and  Ti(0)  =  1. 


— Q  5 


Ti(t)  =  ^6ne'"‘, 
n=0 


Now  consider  the  non  homogeneous  equation 


^  G{t  -  t)  +  HT2{t'^  dr  +  +  ^T2{t)  =  J^{t) 


[  T2(0)  =  0. 


(3.4) 


(3.5) 


By  using  the  transfer  function  and  Lemma  A.l  in  the  Appendix,  we  have 

00  .1 


72(0  =  E°n/  e^'‘('-^)jr(r)dr. 

n=0 


Thus  the  solution  of  equation  (2.14)  is  given  by 


(3.6) 
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(3.7) 


T(f)  =  T,(t)  +  T^it)  =  £  by-'  +  f;  <i„  /' 

„=0  n=0 


Rewrite  ^(t)  as 


J^(t)  =  ^(oo)  -HT  — e-""'  -  2J2  f;  — 6-'^-'-“, 

"  ^2n-l 


n=l 


where 

^■(oo)  =  /f.,!  +  2  J^1)2  +  i  +  iifi  + 

By  (3.7),  we  have 

T'(t)  =  f:  A„6„e^"‘  +  £  g„A„  f  e^''(‘-^).F(r)dT  +  :^(t)  £  a„. 

n=0  n=0  n=0 

Substituting  the  expression  of  T{t)  in  T'(t),  we  obtain 


(3.8) 


where 


nt)  =  E  C'.e^'  +  E  +  E  Cie 

n=0  m=l  m=l 


3  — 1  ^ 


(3,9) 


Cj  =  A„i)„  +  7'(co)a„  -  /fa„A,  E 

m= 


"I”  ^m) 


-2j2a„A„  ^ 


=  1  ^2m— l(An  ”h  ^2m— l) 


Cl  =  ffEr-Tr^TT-i^E-.; 


n=0  +  '^m)  „=o 


d  =  2j^y: 


~  anA„ 


2J2 


G 

n=0  ^2m— l(^n  ^2m— l)  ^2m  — 1  ^=0 


Let  US  first  consider  Simplifying  C^,  we  obtain 


ci=-ht:^ 

„_0  +  c„ 


Let  s  =  — Um  in  (A.l),  we  have 


-Vm  -  A„ 


=  0. 
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Hence  =  0.  By  the  same  reason,  we  also  have  =  0.  Let  us  estimate  C^.  Simplifying 
the  expression  of  and  substituting  (3.8),  we  obtain 


cl  =  XnK  +  +  cxj)an  +  HanG{Xn)  +  2J^a„  X  ,  •  (3.10) 

Z  rn=l  ^2m-l 


By  Lemma  A.2  in  the  Appendix,  substituting  6„  =  (G(A„)  +  /i)an  into  (3.10)  and  using  the 
fact  that  <?(A,„)(A,„  +  H)  +  liX^  + =  0  for  any  m  >  0,  we  obtain 

1  °°  1 

c\  =  ^{pdJ^  +  <^JK  +  2 x;  ,  (3.11) 

I  A„  +  V2m-1 

Thus 

nt)  =  E  E  '  (312) 

n=0  m=l  +  ^2m-l  / 


Let  J  =  —  l^l,  where  |J|  is  the  magnitude  of  J.  Suppose  there  is  to  >  0  such  that 
T'(to)  =  0,  then  it  is  easy  to  check  that  T"(to)  >  0.  Hence  T{t)  may  have  at  most  one  local 
minimum  point  and  no  maximum  point.  By  solving  T'(to)  =  0  for  [  Jj,  we  have 


1J|  = 


oc  ^  ^  n^e 

n=0 


A„to 


(3.13) 


o^nto 


P<iE<-.^'-'-+4E“«  ErrrE 

71=0  71=0  \m=l  ^  '  2m  1 J 

When  to  increases,  \J\  decreases.  Note  that 

•  ••  <  A„  <  A„_i  <  •••  <  Ai  <  Ao  < -if,  and  A„  ~  ^(n^). 


Then  by  letting  to  oo,  we  have 


1.71  = - - .  (3.H) 

Equation  (3.14)  is  a  nonlinear  algebraic  equation  because  Ao  is  dependent  of  J.  The  exis¬ 
tence  of  solution  to  (3.14)  is  easily  obtained  by  using  the  Mean  Value  Theorem,  while  the 
uniqueness  follows  from  the  monotonicity  of  Ao(/)  and  Role’s  Theorem.  Denote  the  solution 
to  (3.14)  by  Jq.  When  J  >  Jq,  equation  (3.13)  will  admit  a  solution  for  to,  hence  T'(to)  =  0. 
Thus  the  solution,  T(t),  to  (2.14)  with  J  >  Jq  decreases  first,  reaches  its  minimum  value  at 
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t  =  to  and  then  increases  to  its  steady  state  Te©.  On  the  other  hand,  the  solution  T{t),  to 
(2.14)  with  J  <  Jq  decreases  monotonically  to  its  steady  state  Too  because  T'(t)  7^  0  for  any 
finite  value  of  t  >  0. 


Next  we  will  give  a  lower  bound  estimate  of  the  value  Jq.  Notice  that 


00  1  00  1 

^ 

,^1  A„  +  V2m-l  Ao  +  V2m-1  ’ 


Vn  >  0. 


Thus  form  (3.11),  we  have 


c;  <  kpdlJI^  -  a|J|)a.  +  2|J'pa„  £  - .  (3.15) 

z  Ao  +  V2m-1 

We  need  to  estimate  the  value  of  Aq.  Let  s  =  —II  —  0(vi  —  H)  where  0  <  0  <  1  to  be 
determined.  Then 


G{s){s  +  H)  +  IJ.S  +  u 

oo  -i  n 

e 

< - -  + —  iiH. 

1-e  ^ 

By  solving  the  following  inequality 

we  obtain  that 

^  1/  —  iiH 

9  > - - - . 

-  I 

Thus  by  letting 

Q_ 

\  +  v-nH' 

we  have 

-n<-H-  -H)<X„<  -H.  (3.16) 

1  +  1/  —  nH 

°°  1  .  . 

Since  ^ -  is  a  continuous  and  decreasing  function  on  (— Vi,0),  it  follows  from 

m=l  ^  +  ^2m-l 

(3.15)  that 
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c‘  <  a„  +21Jp  £  v-ytH  , 


pd\  .|2  <^|  T|  I  2|Jp  ^  1 

=  yi-^i  — ;;  y 

^  ^  m=i - tl — —  +  (2m  -  if 

\  1  +  1/  —  nH 


TT 

/  TT 

+■  cj  n  1  —  i 

\ 

V2\ 

f  1  +  1/  -  pH ) 

/ 

,pd,„<)  a,,,  2|Jp/\/37ry  /' tt  /  u  -  ^H 


where,  in  deriving  the  last  inequality,  we  have  used  the  fact  that  <  1  and  v  —  fiH  > 
By  using  the  inequality 


:,-,iH=\a\J\  +  \  +  ^{l-C:), 


then  we  obtain 


/^l  rl2 


{a\J\^  +  b\J\-c), 


where 


b  =  pd+  —  +  ^[l  +  Hd{l-Ct)], 

TT  TT** 


c  =  a. 
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Let 


T*  — 
Jn  - 


—6  +  +  4ac 


^  2a 

then  <  0,  Vn  >  0  when  0  <  |  Jj  <  Jq.  Hence  from  (3.12),  T'(t)  <  0  for  t  >  0. 


Renicirk  3.3  Corresponding  to  the  example  given  in  the  last  section,  with  the  data  given 
in  Table  1,  we  have 


Jq  «  0.9780. 

Our  computation  (see  Fig.  2)  confirms  that  when  \  J\  <  Jq,  T{t)  is  always  decreasing,  and 
for  some  values  |  J]  >  Jq,  T{t)  decreases  first  and  then  increases  to  its  steady  state  T^o. 


Remark  3.4  The  results  in  Theorem  3.2  cease  to  hold  whenever  is  a  function  of  tem¬ 
perature,  which  is  the  case  during  the  phase  transformation  of  SMA.  Even  though  we  have 
not  proved  such  a  theorem  yet  for  the  variable  case,  a  similar  result  should  be  valid  for 
the  variable  Q  case.  In  fact,  the  constant  case  is  an  approximation  of  the  variable 
case,  hence  the  formula  for  the  value  Jq  may  be  used  in  the  variable  case  approximately 
by  replacing  by  an  average  <C^>.  This  is  numerically  confirmed  by  the  example  given 
in  this  paper,  as  shown  in  Fig.  3,  for  C"{T)  given  by  (2.19)  and  (2.20). 

4  Conclusion 

In  this  paper  we  studied  the  solution  behavior  of  the  transient  heat  transfer  problem  of 
the  1-D  symmetric  thermoelectric  shape  memory  alloy  (SMA)  actuator.  From  the  analytical 
model  proposed  in  [13],  we  proved  that  for  the  transient  cooling  problem  with  constant 
electric  current  density  of  magnitude  |  J|  there  is  a  value  Jq  of  |  J|  such  that  when  j  J]  <  Jq, 
the  temperature  in  SMA  is  always  decreasing  to  its  steady  state,  and  when  |J|  >  Jq,  the  , 
temperature  in  SMA  may  not  be  always  decreasing,  which  is  an  important  property  of 
thermoelectric  SMA  actuators.  A  lower  bound  for  the  value  Jq  was  obtained. 

Due  to  the  relatively  large  heat  conduction  coefficient  of  the  SMA  compared  with  the  heat 
conduction  coefficient  of  the  P-type  and  N-type  semiconductors,  the  temperature  distribution 
in  the  SMA  material  is  expected  to  be  almost  constant  for  small  ratio  d/L,  and  the  results 
derived  in  this  paper  are  valid.  When  the  ratio  d/L  is  not  small  and  the  physical  domain 
of  an  SMA  actuator  is  no  longer  1-D,  a  proper  three  dimensional  model  of  SMA  actuators 
should  be  introduced  and  studied. 
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Appendix 

Lemma  A.l  Let  J{t)  =  J  <  0.  Let  w(s)  be  the  transfer  function  of  equation  (2.14).  Then 
there  exists  a  sequence  of  positive  numbers  {a„}^o 


i&(5)  =  -t: - =  5  G  3?, 

(?(s)(5  +  H)+  ps  +  u  -  Xn 


where  {A„}^o  single  poles  ofw{s). 


(A.1) 


Proof.  From  the  proof  of  Theorem  3.1,  we  have 


•  •  •  <  A„  <  A„_i  <  •  •  •  <  Ai  <  Ao  <  -if,  and  A„  ~  O(n^), 


because  J  <  0.  Generalize  the  transfer  function  w(s)  to  the  complex  plane  C,  i.e.  w{z), 
z  G  C.  Thus  w{z)  is  analytic  everywhere  except  {An}  and  each  An  is  a  simple  pole  of  w{z). 
The  residue  of  w(z)  at  z  =  A„  is  given  by 


Res{w{z),  A„) 


r  00 


E 


Vm-H 

{Xn+Vm? 


1-1 


+  P 


Let  a„  =  Res(w(z),Xn).  Then  On 
because  A„  c;i  O(n^).  Let 


>  0  and 


OO 


E 


z  -  A„ 


is  convergent  for  any  z  €  C  \  {A„} 


f(z)  =  w(z)  - 

n=0 

Since  /(z)  has  no  poles  in  C  and  ih(^)  is  analytic  on  C  \  {A„},  /(z)  is  an  entire  function  [7]. 
We  will  prove  that  /(z)  =  0,  z  G  C. 

Let  z  =  {vm  —  H)e'^  —  H,  6  £  [0, 27r]  and  m  be  sufficiently  large.  We  have 


lu;(z)| 


_ 1 _ 

\G{z){z  +  H)  +  pz  +  u 
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me 


HT^m 


+  n^  e’^  +  1 


< 


e.-5  .  1  +  -{^-f^H)-J2 

*  n^m 


|n^  —  m^l 


By  the  inequality  |a  +  6i|  >  |a|  for  any  a,  5  6  3?,  we  obtain 


1 

I"  ^ 


e‘®  +  l  C„I2 


1  9  0 

>  -  +  fim  Tt  . 


Note  that 


V  ^ 

n^m  1^^  - 


m-1  ^2 


dx  + 


m" 
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m" 

+  /  - 7:dx 


2  ~  (m  ~  1)2  (m  +  1)2  —  m2  Jm^i  x^  —  m2 
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4m^  m  ,  m  +  X 


+  —  In 


2m^  —  1  2  m  —  X 


m— 1 


m  ,  X  —  m 


+  —  In 
0  2  X  +  m 


m+l 


4m^  m  ,  , ,  m  ,  _  , , 

+  _  ln(2m  -  1)  -  -  ln(2m  +  1) 


<  2m +  2. 
We  hence  obtain 


\H^)\  <  Y 


+  jj.m^T:'^  —  {u  —  IJ.H)  —  (2m  +  2) 


Also  for  z  =  (vm  —  H)e'^  —  H,  6  £  [0,27r]  and  m  >  2,  by  using  the  properties  of  {A„}  and 
the  expressions  of  {a„},  we  have 


OO 

n=0 


<  E 


folium -H)e'^-iX,+H)\ 
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<  E 


^0  l^m  +  A„| 

m— 1 


<  —  1  _| _ _ 
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1  /2m  +  2\  4 
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\  J  m 


Therefore  for  z  =  {vm  -  H)e'^  —  H,  6  e[0, 27r]  and  m  be  sufficiently  large,  we  have 


|/(z)|  <  |u)(2)|  + 


^  T 


V' 

Z  —  Xn 
1 


1  2Tn  +  2  4 

4 - 1 - . 

2  4112  yjj 


-  +  /imV^  -  (i/  -  -  (2m  +  2) 


Since  /(z)  is  an  entire  function,  by  applying  the  Maximum  Modulus  Theorem  [7],  we  obtain 
that  /(z)  =  0  on  C.  Hence  (A.l)  is  proved.  I 

By  the  same  way,  the  following  lemma  can  also  be  derived. 

Lemma  A.2  Let  J{t)  —  J  <  0.  Then  there  exists  a  sequence  of  real  numbers  {fcn}^o  such 
that 

_ M  +  ^(g) _ _  bn  5  £ 

G(s)(s  +  H)  +  JJlS  +  u  „=o  ^  An 

Furthermore, 

OO 

bn  —  (G(An)  4"  ff)an,  and  'X  .  bn  1. 

n=0 


(A.2) 

(A.3) 
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